Gravitational collapse without a remnant 
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We investigate the gravitational collapse of a spherically symmetric, inhomogeneous star, which 
is described by a perfect fluid with heat flow and satisfies the equation of state p = p/3 or p = Cp 1 
at its center. Different from the ordinary process of gravitational collapsing, the energy of the whole 
star is emitted into space. And the remaining spacetime is a Minkowski one at the end of the 
process. 

PACS numbers: 97.60.-s, 04.40.Dg, 04.25.Dm 



I. INTRODUCTION 

Gravitational collapse is one of the most impor- 
tant topics in general relativity. A large mount of 
models of gravitational collapse has been proposed 
[1-7]. In these models, the remnant of a gravita- 
tional collapse is possible a compact self-sustained 
star or a black hole, or a naked-singularity, de- 
pending on the initial condition of the collapse. A 
decade ago, in studying the critical phenomenon 
in gravitational collapse, first revealed by Chop- 
tuik [5], Hirschmann and Eardley showed that 
the gravitational collapse of a complex scalar field 
may leave behind an approximately flat space- 
time. Their model about the collapse, as a so- 
lution of Einstein field equations, possesses self- 
similarity [6]. Later, Schafer and Goenner con- 
structed a model of gravitational contraction of a 
radiating spherically symmetric body with heat 
flow [7], in which both initial mass and initial 
radius are infinitely large and all mass of the 
body can be radiated away eventually without 
forming an event horizon. Also about a decade 
ago, Fayos, Senovilla and Torres studied geometry 
matched by two spherically symmetric spacetimes 
through a timclike hypersurface from a very gen- 
eral point of view and exhausted all possible and 
qualitatively different matchings with their corre- 
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sponding conformal diagrams for a flat Robertson- 
Walker model with a linear equation of state p = 
7/0 [8]. In particular, Fig. 9 and the time rever- 
sal of Fig. 11 in [8] show that a flat spacetime is 
left when a radiating homogeneous star (or Uni- 
verse) with a linear equation of state radiates its 
all mass. 

In the present paper, we propose a new ap- 
proach to study the gravitational collapse of a 
spherical star. In our approach, a fluid star with 
finite initial mass and radius is supposed to be 
spherically symmetric, inhomogeneous, with the 
equation of state p = p/3 or p = Cp 7 at its cen- 
ter, and having heat flow outside it. Similar to 
the solutions given by Hirschmann and Eardley [6], 
Schafer and Gocnner[7] and Fayos et al[8], all en- 
ergy of the star in our model will be emitted in 
the process of collapse, and the remaining space- 
time is an empty flat one. It is remarkable that 
for a star with about a solar mass and a solar 
radius, the energy at the order of 10 54 erg will 
be emitted into space within about 8 ~ 600s. 
As a result our model concludes that the inferred 
isotropic average luminosities are on the order of 
10 51 ~ 53 erg/s, which has the same order of magni- 
tude for a gamma-ray burst. 

The arrangement of the paper is as follows. The 
method of model construction and the equations 
governing the model arc described in next section. 
The junction conditions and boundary conditions 
are listed in Section 3. In Section 4, as examples, 
several numerical solutions are presented. The 
concluding remarks are given in Section 5. 



II. FORMULATION 

In the pioneer paper to understand the late 
stages of stellar evolutions [1], a collapsing star 
is supposed to consist of homogeneous, spheri- 
cally symmetric, pressure-less, perfect fluid and 
to be surrounded by an empty space. The interior 
of the star may be described by the Friedmann- 
Robertson- Walker metric [9] 

dr 2 



ds z 



dt 2 - a 2 {t) 



r l dVL' 



(1) 



1 — kr 2 

where d£l 2 = dO 2 + sin 2 8dip 2 is the metric on a 
unit 2-sphere, k is pp times of the initial energy 
density of dust in the unit of c = 1, and a(t) is the 
solution of 

a 2 (t) = k[a- 1 (t)-l\. 

In an astrophysical environment, a star usually 
emits radiation and throws out particles in the 
process of gravitational collapse. In this situation, 
the heat flow in the interior of a star should not 
be ignored and the exterior spacctime is no longer 
described by a Schwarzschild metric. To take the 
radiation of a star into account, the interior solu- 
tion of the gravitational collapse of radiating stars 
should match to the exterior spacetime described 
by the Vaidya solution [10] 



ds 2 = (1 - 



2GM(v) 
R 



)dv 2 + 2dvdR- R 2 dn 2 , (2) 



which has been studied extensively [2-4, 8]. In 
particular, the gravitational collapse of a radiating 
spherical star with heat flow has been studied in 
an isotropic coordinate system [4] 



ds 2 =dt 2 - B 2 {t,r)(dr 



? -dn 2 ), 



(3) 



where 
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with suitable b(t) and A(t) functions. They called 
the solution the Friedmann-like solution. 

We study the gravitational collapse of a spheri- 
cally symmetric, inhomogeneous star in a proper- 
time reference, which can be written in a general- 
ized Friedmann coordinate system 

dr 2 



ds' 



dt 2 -a 2 {t) 



1 - k{t)r 2 



r 2 dtt 2 



(4) 



with k being a function of t. 

The stress-energy tensor of the fluid with heat 
flow and without viscosity is given by 

T^ = {p + p)u^u v - pg^ + g"u" + q v u ti , (5) 

where p and p are the proper energy density and 
pressure measured by the comoving observers re- 
spectively, u^ is the 4- velocity of the fluid, and g M 
is the heat flow. u M and g M satisfy 

u M u" = 1 (6) 

«„u" = 0. (7) 

For the spherically symmetric collapse, one has 

u" = (uV\0,0) (8) 

cf = (gV,0,0). (9) 

The Einstein's field equations 

G^ = -8nGT^ (10) 

and the covariant conservation of stress-energy 
tensor give rise to 

o r, k + d 2 + 2aa Y 

8-n-Gp = 3 + - ir - r >, (11) 



lirGp 



2 ' 2 a2 ' 

a a A 

k + a + 2ad X 



a 



a 2 A 2 



, (12) 



"(f 



«i 



\ 
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8irG qi = 1 
with the identity 
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(14) 

(15) 
(16) 
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where the over-dots denote the derivatives with 
respect to time t, and 

1 



y/l-k(t)r 2 
X = 3aaAA + a 2 AA, 



(18) 
(19) 



Y [2(k + d 2 -atl)A-aciA-a 2 A]A, (20) 

(21) 
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X, Y, Z should satisfy the following solvable con- 
dition, 



Z l > 4XY. 



(22) 



For spherical collapse, the dominant energy con- 
ditions [11] require 



P-P>0, 

p + p>2^/-q^, 



which lead to 



X + Y > \Z\ >0. 



III. BOUNDARY CONDITIONS 



(23) 
(24) 



(25) 



At the center of a star r 
thus X, Z vanish. Hence, 



0, A = A = and 



MoIr=0 = lj Ul|r=0 = 0, 

g o |r=o = 0, g 1 | r=0 = 0, 

inGp\ r=0 = 3a _2 (fc + d 2 ), 



87rGp| r= o = —a 2 (fc 



2ad). 



(26) 

(27) 

(28) 
(29) 



The exterior solution outside the collapsing star 
is described by Vaidya metric (2). The interior 
solution and exterior solution should satisfy the 
Darmois junction conditions on the surface of the 
star [12, 13], 



R s = 


ar s , 


M(v) 


a 
= 2G 


vl = 


q^ l \ s 



(k + a 2 )r c 



(30) 
(31) 

(32) 



and 



[dt 2 - a 2 A 2 dr 2 ] 
2GM{v) 



R 



dv 2 + 2dvdR 



, (33) 



where the subscript s denotes the value taken on 
the surface of the star. In particular, in Eq.(32) 
with the expressions (12), (17), (18), r should be 
replaced by r s (t), for Eq.(32) is only valid on the 
surface of the star. M(v) is the Misner-Sharp 



mass [14] of the star. (Appendix A gives the de- 
tailed discussion on the junction condition.) The 
equation of motion of the surface is 



dr s 
~dt 



- 2 A- 2 ^ 

Uq 



(34) 



In brief, with the help of the equation of state of 
fluid at the center of the star, Eqs.(28), (29), (32) 
and (34) constitute a complete system of equa- 
tions for k{t), ait), r s it). Once kit), ait) and r s (i) 
are obtained, one may directly get the expression 
for pir,t) and p(r, t) from (11) and (12) and then 
determine the Misner-Sharp mass M(v) and the 
surface radius R s from Eqs.(31) and (30). 

IV. NUMERICAL SOLUTIONS 

Without loss of generality, we choose a(0) = 1. 
In order to make numerical calculation conve- 
niently, we introduce the dimensionless physical 
quantities as follows. 



fe- 

A 



kR , 



aRo, 



aR. 



0; 
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AR , A^ARl, p^ p RI/8ttG, 



p -> P Rl/8irG, 



r 



qPRZ/SirG, 



where Ro is the initial radius of star. Note that 
we have already chosen the unit of c = 1. 

In the following, we study several examples nu- 
merically 

A. p = p/3 at center 

Suppose that the equation of state at the center 
of a star take the radiation form, i.e. p\ r =a = 
p/3\ r =o, which was also discussed in Ref. [16]. We 
are interested in the case that the surface of a star 
is almost stationary in the initial state, namely, 

r s (t = 0) = R s it = 0) = R , (35) 

fl,(t = 0)«0. (36) 

Therefore, we consider the initial conditions (35) 
and 



d(< = 0) = 0, 



(37) 
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FIG. 1: Numerical solution of the gravitational col- 
lapse of spherical star with heat flow in the case 
p| r=0 = p/3| r= o and the Uq solution in Eq.(I3). The 
horizontal axis is the time t in the unit of Ro/c. In the 
process of evolution, r s (t) (dash curve, in the unit of 
Ro) decreases monotonically while a(t) (dotted curve) 
keeps almost a constant. 2M/R„ (real curve) in the 
unit of c 2 /G increases first and then goes to as 
R s — > 0, which implies that the star disappears at 
the end of the collapse. 



FIG. 2: Evolution of p and p at the boundary. The 
horizontal axis is the time in the unit of Ro/c. p and 
p are in the unit of c 4 /(8tyGRq). 



Furthermore, our numerical analysis shows that 
the star with p\ r =o = op| r =o will radiate its whole 
mass in the process of the collapse, without the ap- 
pearance of a horizon, for different initial values. 



In such a case, only the Uq solution in Eq.(13) 
can reach (36) from the initial conditions (35) and 
(37). Fig. 1 presents the numerical solution in the 
case. Since p and p are determined by Eqs.(ll) 
and (12), the initial state of the star is not ho- 
mogeneous and the equation of state in the star 
is generally deviated from the radiation. In the 
figure, the time and radius r s are in the units of 
Ro/c and Rq, respectively. In the process of evo- 
lution, r s (t) (dash curve) decreases monotonically 
while a(t) (dotted curve) almost keeps a constant. 
2M/R B (real curve) in the unit of c 2 /G increases 
first and then goes to more quickly than R s itself 
as R s — > 0. It implies that the star will disappear 
without the appearance of a horizon in the pro- 
cess. 

Fig. 2 shows the evolution of p and p at the 
boundary, in unit of c 4 /(8ttGRq). The emission 
of the star arrives its maximum value at about 
0.9Rq/c after the beginning of the collapse. In 
the late stage of the process, the equation of state 
at the boundary as well as at each point in the 
star tends to p = \p. At the end of the process 
both the energy density and the pressure become 
0, which confirms that the whole star is radiated 
out into space in the process. 



B. 



= Cp 1 at center 



Now, let us suppose that the equation of 
state at the center have the forms of poly- 
tropes. Again, we only consider the u^ solu- 
tion in Eq.(13) with the initial conditions (35) 
and (37). Figs. 3 and 4 present the the numer- 
ical solutions of Misncr-Sharp mass in the cases 
7 = 6/5, 11/9, 5/4, 9/7, 4/3, 7/5, 3/2, 5/3 and 
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FIG. 3: The ratios of mass to radius of collapsing 
stars for numerical solutions with p| r =o = p 7 |r=o, and 
7 = 9/7, 5/4, 6/5, 11/9. t is in the unit of Ro/c. 
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FIG. 4: The ratios of mass to radius of collapsing 
stars for numerical solutions with p\ r =o = p 7 |r=o, and 
7 = 5/3, 3/2, 4/3, 7/5. t is in the unit of Ro/c. 



FIG. 6: Evolution of p at the boundary for the cases 
of 7 = 5/3, 3/2, 4/3, and 7/5. 



C = 1 , respectively. The time is again in the unit 
of Ro/c. The behaviors of r s for all cases are sim- 
ilar to the one in Fig.l. These figures show that 
the evolutions of the stars with different adiabatic 
indexes have the similar behavior to that in Fig. 1 
but with different collapsing time scale. Namely, 
the whole star will be emitted out without a rem- 
nant eventually in all these cases. 

Figs. 5 and 6 show the evolutions of p at the 
boundary for p — p 1 with 7 = 9/7, 5/4, 6/5, 
11/9, 5/3, 3/2, 4/3, and 7/5, respectively, in unit 
of c 4 /(8ttGRq). Figs. 7 and 8 show the evolutions 
of p at the boundary for the same parameters. 



The times to emit all energy, T, are shown in the 
following Table. The last line in the Table gives 
the numerical values for the star with a solar 
radius initially. At the end of the process both 

Table. Time to emit all energy for the star with 
p = p 1 at its center. 
7 6/5 11/9 5/4 9/7 4/3 7/5 3/2 5/3~ 
cT/Rq 11 14 18 27 45 94 262 630 
T(s) 25.3 32.7 42.0 62.7 104.5 218.7 608.2 1462.5 

the energy density and the pressure become 0, 
which confirms that the whole star is radiated out 
into space in the process. 

Finally, we give the evolution of k(t) and a in 
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FIG. 5: Evolution of p at the boundary for the cases of 
7 = 9/7, 5/4, 6/5, and 11/9. The horizontal axis is the 
time in unit of Ro/c. p is in the unit of c 4 /(8tvGRq). 



FIG. 7: Evolution of p at the boundary for the cases of 
7 = 9/7, 5/4, 6/5, and 11/9. The horizontal axis is the 
time in unit of Ro/c. p is in the unit of c 4 /(8ttGRo). 
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FIG. 8: Evolution of p at the boundary for the cases 
of 7 = 5/3, 3/2, 4/3, and 7/5. 



k(t)*10 b 




10 12 



FIG. 9: Numerical solution for k(t) and a(t) in the 
case of p\ r =o = P |r=o and the Uq solution of 
Eq.(13). 



Fig. 9. Obviously, k(t) and a have the same order 
of magnitudes, which is consistent with the fact 
that the equation of state at the center of the star 
serves as a free input parameter. 

V. CONCLUDING REMARKS 

We have shown a new approach to study the 
gravitational collapse of a spherical, inhomoge- 
neous, fluid star with heat flow. By use of the 
Ansatz of the generalized Fricdmann-Robertson- 
Walker metric, we obtain the new solutions of 
Einstein's field equations. These solutions share 
the same properties: the initial mass and radius 



are finite and whole stars will be emitted in the 
process of the gravitational collapse in a finite in- 
terval, without the formation of a horizon, and 
a Minkowski spacetime is left at the end of the 
process. This is quite different from the standard 
evolution of the gravitational collapse of a star, fn 
the standard process, a white dwarf, or a neutron 
star, or a black hole, or even a singularity will be 
formed at the final state and only a part of whole 
mass will be emitted. Obviously, the numerical so- 
lutions of the Einstein's field equations presented 
here are the description of another type of collaps- 
ing evolution beyond of the standard process. 

As mentioned in Introduction, the numerical so- 
lutions with the property that all energy of the 
star in our model will be emitted in the process of 
collapse and an empty flat spacetime will be left 
behind have been obtained before [6-8] . Neverthe- 
less, our solutions arc more realistic than the pre- 
vious ones. This is partly because in our solutions 
there is no self-similarity which will break down 
the asymptotic flatness of spacetime and partly 
because there is no initial singularity and the ini- 
tial mass and radius are finite. Another reason is 
that the interiors of stars in nature are inhomo- 
geneous and having the equation of state of poly- 
trope. Our solutions presented here share these 
characters. In particular, when we apply our so- 
lution to a star with about a solar mass M Q , a 
solar radius Rq and p = p/3 at the center, a huge 
mount of energy (about 1.8 x 10 54 erg) will be 
emitted into space within 7.57 seconds, which are 
typical values for a gamma-ray burst. [17] If stars 
have the same size and p — p 1 at their centers, 
the times to emit all energy, T, will range from 
25s to 1462.5s. They are all reasonable values for 
gamma-ray bursts. And according to our numeri- 
cal solutions we believe that if 7 changes continu- 
ously, then the times to emit all energy will change 
continuously. Thus, the gravitational collapse of 
such a star might provide a new energy mecha- 
nism for gamma-ray bursts. Of course, detailed 
investigations of applying the solution to describe 
gamma-ray bursts are needed. 
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APPENDIX A: BOUNDARY CONDITIONS 
ON THE STAR SURFACE 



The motion of the star surface can be described 
by a time-like three-space E which divides the 
spacetime into the interior V~ and the exterior 
V + . The induced metric on and the extrinsic cur- 
vature of E are 



d 4| E 



{g^ dx ± dx ±) 



(Al) 



K± = -r£ 



, , dx+ dx+ , . , 
n±r±?- --^-, (A2) 



respectively, where x± , g^ v and T^" are the coor- 
dinates, metrics and affine connections of V re- 
spectively; t; 1 is the intrinsic coordinate of S and 
n„ are the unit covariant vectors normal to E. 

Following Israel [15] , the junction conditions are 
the continuity of the induced metric and extrinsic 
curvature, 



rfs -ls= ds +| S ' 



Kl 



Kt 



(A3) 
(A4) 



For the spherically symmetric collapse with 
heat flow and a radiating surface, a detailed dis- 
cussion was given by Santos [13] with the interior 
solution in isotropic, comoving coordinates 



dst = A 2 (t, r)dP - B 2 (t, f)[df 

and the exterior Vaidya solution 
2 _ n 2GM(v) 



*Wl 



(A5) 



ds% = (1 



R 



)dv 2 + 2dvdR - R 2 dn 2 . (A6) 



A similar discussion to Santos' can be taken for 
the more general interior solution in comoving co- 
ordinates, 

dsl = e^dP - e^dr 2 - TZ 2 dn 2 , (A7) 



where <p, ip and TZ are functions of i and f . 

In the comoving coordinates, the equation 
of motion of the star surface is f = fs(= 
constant). The corresponding hypersurface equa- 
tion is f~(i, f) = r — fs = 0. The space-like vec- 
tor df~/dxt is normal to E and by normalization 
g- U n~n~ = — 1 one will have 



n- = (0,e*,0,0) 



(A8) 



By utilizing the metrics (A6) and (A7) in the first 
junction condition (A3), one can have 



R\^ = n(t,fs), 



(A9) 



[(1 - 2GM/R)dv 2 +2dvdR]^ = e^dt 2 ^ .(A10) 

The equation of motion of star surface in the exte- 
rior space is given by the above equations. By sup- 
posing the hypersurface equation as f + (v,R) = 
R — R(v) = 0, the normal covariant vector n+ can 
be written as 



av 



,1,0,0). 



(All) 



The normalization g^n^n^ = — 1 gives 



A 2 [(1 - 2GM/R) + 2{dR/dv)] Yi = 1. (A12) 

Utilizing the equation (A10), one easily gets A = 
e~^v\s, and then 

n+ = e _ *v(-ditydu,l,0,0)| s 

= e-*(-A,r} J 0,0)| El (A13) 

where over-dots denotes the derivative with re- 
spect to t. 

The intrinsic coordinates of E can be conve- 
niently chosen as £* = (i,8,ip). Now with the 
obtained normal vectors and metrics, the induced 
extrinsic curvature can be calculated out. Their 
nonzero components are 



K J 

Kl 



K 






K$ 



K+ 



Kl a sin 2 6 

V 

V 



'R 



R 



GM 



n 2GM v 



= K ee sin 



(A14) 
(A15) 
(A16) 

(A17) 

(A18) 
(A19) 



The second junction condition (A4) (having only 
two independent equations because of spherical 
symmetry) can therefore be rewritten as 



-e-V| s = e~* 



e^RR'l 



e-VR 



GM 



■ , 2GM S 
R+(l-—-)v 



(A20) 

(A21) 



With the help of the metric continuity equations 
(A9) and (A10), the equation (A21) can be re- 
duced to 



M{v) 



R 
2G 

m(t, r) 



1 + e -20^2 _ e -20^/2 



(A22) 



where m(t, f) is the well-known Misner-Sharp 
mass function [14]. Substituting M(v) into the 
equation (A21), one also gets 



* 



e-*TZ + e-^W 



(A23) 



Substituting M(v) and v into equation (A20), one 
now gets 

e-t-^CRfi + R'ip-R') 



e -24> (ji + K_^ 



2R 



-74 



'RP 

2R 



K '*'' + 2K 



(A24) 



Suppose that the collapse material be a shear- 
free fluid with heat flow described by the stress- 
energy tensor (5). In the comoving coordinates, 
the unit four-velocity vector u^ = (e~*, 0,0,0) is 
normal to the space-like heat flow vector g M = 
(0, g 1 , 0, 0). By solving the Einstein equation, one 
gets 



BTrGg 1 
8nGp 



-<p-2ip 



R 



(R,4>' +R'ijj-R,'), (A25) 



,-20 



■■ R 2 
— \R(j)-R-—, 



+e 



R 



-2tl- 



R 2 



'R' 2 2 
R 2+ R 



R' 



R 2 



(A26) 



Now from the equations (A24), (A25) and (A26) 
it is easy to see that 



Pis = 



<t q» 



(A27) 



In the end, one can conclude that for spheri- 
cally symmetric collapse with heat flow but shear 
free, the boundary conditions are equations (Al), 
(A22) and (A27), which are all invariant for the 
coordinate transformations on time-radial plane. 
Therefore although the coordinate is not comoving 
in our model, the above three boundary equations 
are still used directly. 
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